
NASA 

Technical 

Paper 

2452 

1985 


Slip-Boundary Equations 
for Multicomponent 
Nonequilibrium Airflow 


Roop N. Gupta 

Old Dominion University 
Norfolk , Virginia 

Carl D. Scott 

Lyndon B. Johnson Space Center 

Houston , Texas 

James N. Moss 

Langley Research Center 
Hampton , Virginia 


l\J/\SA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 



Acknowledgment 

The authors are grateful to Dr. Fred G. Blottner, of Sandia National Laboratories, 
Albuquerque, New Mexico, for many useful discussions pertaining to the mass flux 
expressions contained in appendix B. 



Contents 

Symbols v 

Summary , . . , 1 

Introduction . . 1 

Analysis 1 

Various Fluxes and Distribution Functions , . . . 3 

Balance and Slip Equations 5 

Simplifications for Multicomponent Mixture 11 

Concentration Slip Boundary Condition for Fully Catalytic 
or Noncatalytic Surface .............................. 14 

Simplifications for Binary Mixture 14 

Comparison With Existing Expressions 16 

Simplifications for Single-Species Gas 18 

No-Slip Boundary Condition for Species Concentration 19 

Discussion . . 21 

Concluding Remarks . . .21 

Appendix A — Constants Appearing in Distribution Function and 
Diffusion Velocity in Terms of Transport Properties . . . 22 

Appendix B — Evaluated Expressions for Fluxes of Mass, Momentum, 

and Energy From Integrals Over Velocity Distribution Function 24 

Appendix C — Slip Boundary Conditions Given by Davis for Binary Mixture . 29 

Appendix D — Slip Conditions in Body-Oriented Coordinate System 31 

Appendix E — Slip Conditions in Spherical Coordinate System 42 

Appendix F — Integrals Required for Obtaining Various Fluxes 46 

References 50 

Figures 

1. Knudsen layer showing general fluxes and coordinate axes; 

temperature as function of normal distance is schematically overlayed . 2 

Dl. Coordinate configuration for body that is symmetric about its axis 31 

El. Spherical coordinate configuration 42 




Symbols 


4- He 4 W 

(j) 

a iQ- a i\i b{Qi c j'q 


b 0 
Ci 

Cp 


^p,t 

Cv 

D ij 

Of 

Dam 


J 

a k 

E 

E t 


e 


F 

f D 

r 

f w 


h 

hi 


A hi 
K 
k 

Ir 

Wl 


,298 


coefficients of velocity distribution function for *th species 


distribution function coefficients given in appendix A and 
defined in reference 7 for ith species 

speed of sound, \f^RToo 

distribution function coefficients defined by equation (70c) or 
(DIO) 


mass fraction of species i 

frozen specific heat at constant pressure of mixture, 



specific heat at constant pressure of species i 

frozen specific heat at constant volume 

multicomponent diffusion coefficient for species pair i and j 


thermal diffusion coefficient 

binary diffusion coefficient for mixture of atoms and 
molecules 

binary diffusion coefficient for mixture of species i and j 


diffusion vector of jth species 
total energy flux, E ^ + eM 
energy flux from translational energy 

internal energy that readily equilibrates with translational 
energy 

general flux of property such as mass, momentum, or energy 

diffusion coefficient, P^Um/M^oo/^oo) 

distribution function of ith species 

equilibrium (or Maxwellian) distribution function of zth 
species 

enthalpy of mixture, Cjh^ 

T , 

enthalpy of species i, f c v MT + Ah - 9Q8 
298 ’ 

heat of formation for species i at 298° K 


thermal conductivity of mixture 
Boltzmann constant 

wall catalytic recombination rate constant for species i 


v 



M 

mass flux 

Moo 

Mach number, U 0 0/0,00 

m i 

mass of ith. species 

w, 

mass of mixture or single-species gas molecule, 1 ~ j 

N Le 

binary Lewis mrnaber, CpP$ AM ! K 

^Le,ij 

multicomponent Lewis mnnber, CppD^j/K 

^ Le,ij 

binary Lewis number, CppO^j/K 

N Pr 

Prandtl number, Cpp/K 

NS 

number of chemical species 

n i 

number density of species i 

n* 

coordinate measured normal to body in body-filled coordinate 
system (appendixes C and D) 

ft, 

total number density, 

P 

momentum flux 

P 

pressure 

Pi 

partial pressure of species i 

Qw 

wall heat-transfer rate or wall heat flux 

R 

universal gas constant divided by molecular weight of mix- 
ture, R/W 

R 

universal gas constant 

R* 

radius measured from axis of symmetry to point on body 
surface 

r 

radial distance from body center (appendix E) 

r N 

body nose radius of curvature (appendix D) or radius of 
sphere (appendix E) 

s* 

distance measured along body surface (appendix D) 

T 

temperature 

Uoo 

free-stream velocity 

u 

velocity component along body surface (appendixes C and D) 

V 

thermal velocity vector 

n 

thermal (or peculiar) velocity of *th species (component of V 
for species i) 

v{ 

diffusion velocity 

4 

total (or molecular) velocity, vq^ + Vfe 

m 

bulk or mass-averaged velocity 


vi 



n 

Wi 

W 

w 0k 

x 

y 

z 

a 

P 

7 

7 

hi 

£ 

6 

K 

X 

M 

P 

$*(V) 

<W) 


*l>i 

Subscripts: 

A 


e 

M 

M 

q 

ref 


x,y,z 


y 


II , x, z 


nondimensional thermal velocity, 



molecular (or atomic) weight of species i 


molecular weight of mixture 


nondimensional mass-averaged velocity v ok 


coordinate parallel to body (fig. 1) 
coordinate perpendicular to body (fig. 1) 
coordinate parallel to body (fig. 1) 
energy accommodation coefficient 
body angle defined in figure D1 

recombination coefficient (or, fraction of incident atoms that 
recombine at surface) 


ratio of frozen specific heats, c p /c v 
Kronecker delta 


Reynolds number parameter, (// re f /poo £Ax> r/v ) 1 / 2 

accommodation coefficient (or fraction of incident particles 
that stick) 

surface curvature 

mean free path 

viscosity of mixture 

density 


perturbation part of distribution function 
circumferential angle (appendixes D and E) 

general property of *th species such as mass, momentum, and 
energy 

source term defined by equation (38) 


atom 

equilibrium value 
coordinate indices 
molecule 
species index 
reference quantity 
component directions 
normal component 
tangential components 


vii 



00 


free-stream value 


Superscripts: 


I 

incident flux 

T 

specularly reflected flux 

Affixes (used both 

as subscripts and superscripts): 

hi 

species indices 

s 

edge of Knudsen layer 

w 

wall 


A bar over a symbol indicates a dimensionless quantity except for the definitions of and V\. 


viii 



Summary 

Equations are presented for the surface-slip (or jump) values of species concentration, pressure, 
velocity, and temperature in the low Reynolds number, high-altitude flight regime of a space vehicle. 
These are obtained from closed-form solutions of the mass, momentum, and energy flux equations by 
using the Chapman-Enskog velocity distribution function. This function represents a solution of the 
Boltzmann equation in the Navier-Stokes approximation. The analysis, obtained for nonequilibrium 
multicomponent airflow, includes the finite-rate surface catalytic recombination and changes in the 
internal energy during reflection from the surface. Expressions for the various slip quantities have 
been obtained in a form which can readily be employed in flow-field computations. A consistent set of 
equations is provided for multicomponent and binary mixtures and single-species gas. An expression 
is also provided for the finite-rate species-concentration boundary condition for a multicomponent 
mixture in the absence of slip. 

Introduction 

For an accurate prediction of the aerothermal environment of a space vehicle entering the Earth’s 
atmosphere in the low Reynolds number, high-altitude flight regime (ref. 1), the multicomponent, 
nonequilibrium gas chemistry, as well as the wall slip and catalysis effects, must be evaluated. Such 
flow fields are of particular interest for aeroassisted orbital transfer and Space Shuttle vehicles. 
This study was undertaken to obtain the boundary relations that incorporate the effects of slip, 
multicomponent diffusion, wall catalytidty, and changes in internal energy of the molecules (during 
reflection from the surface) for application to flow-field calculations under the general assumption of 
local thermodynamic equilibrium. 

Scott (ref. 2) first presented the wall boundary conditions for a multicomponent mixture with 
diffusion and wall-catalyzed atom recombination. In obtaining these boundary conditions, he used a 
first-order velocity distribution function at the edge of the Knudsen layer next to the wall (see fig. 1), 
where the continuum model of the gas breaks down. These boundary conditions, obtained from the 
kinetic theory considerations, provide solutions at the top of the Knudsen layer that would match 
the solution of the Navier-Stokes equations in the bulk outer flow. Hendricks (ref. 3), using Scott’s 
formulation, obtained simplified expressions for engineering applications with some corrections to 
Scott’s expressions. Hendricks’ analysis, however, contained some errors in obtaining the engineering 
expressions. This paper reanalyzes the wall boundary equations by using the approach of reference 2 
and provides appropriate relations for the various quantities with surface slip in a form which can 
readily be employed for flow-field computations. An effort has also been made to reconcile the 
differences between slip expressions employed by the different researchers. The present analysis 
provides a consistent formulation for the slip equations for multicomponent and binary mixtures and 
single-species gas. An expression is also provided for the finite-rate species-concentration boundary 
condition for a multicomponent mixture (in addition to that for a binary mixture) in the absence 
of slip. This expression could readily be used in Shuttle-flow-field heating analyses, which include 
effects of multicomponent diffusion such as reference 4. 

The main difference (among other details) between the results contained in reference 5 and the 
present work is that, in reference 5, the various internal degrees of freedom for a molecule were 
considered frozen during reflection from the surface whereas they were allowed to change in the 
present analysis. 1 Further, the results are provided in both spherical and body-oriented coordinate 
systems. 

Analysis 

The slip conditions are taken to exist across the Knudsen layer, which is on the order of one mean 
free path in thickness as sketched in figure 1. The analysis outlined herein follows the approach of 

1 Since the assumption of local thermodynamic equilibrium is employed, only those internal energies are considered 
which equilibrate readily with the translational energy. 
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F f incident flux; j Fj specularly reflected flux; F™ diffusely reflected flux. 

Figure 1. Knudsen layer showing general fluxes and coordinate axes; temperature as function of normal distance is 
schematically overlayed. 

references 2 and 5. It is based on Shidlovskiy’s (ref. 6) assumption that the distribution function 
near the wall can be described to first-order accuracy by the so-called Navier-Stokes approximation. 
However, a deviation is made from the procedures of Shidlovskiy in that a Chapman-Enskog 
type distribution function for a multicomponent mixture obtained by the variational method of 
Hirschfelder, Curtiss, and Bird (ref. 7) is used. The Chapman-Enskog distribution function permits 
an accurate accounting of diffusion. The analysis contains the following assumptions; 

1. The energy and momentum accommodation coefficients (i.e., a and 6, respectively) have the same 
value. 

2. The fluxes of mass, momentum, and energy across the Knudsen layer are assumed constant. This 
assumption is consistent with the one of negligible variation of the velocity distribution function 
through the Knudsen layer. 

3. The internal energy associated with the rotational and vibrational modes readily equilibrates with 
translational energy. 

The interaction model at the gas-solid interface, with the various fluxes sketched in figure 1, can 
be mathematically stated for dissociated air as provided by the following equations. 2 

2 From the available experimental information, there is no indication that there is significant formation of NO on the 
Shuttle or other thermal protection surfaces. For this reason, the formation of NO has been neglected in the present 
model. 
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For recombining atom: 


Fa = Fa + (» " >a>f\ + (« A - 1a)F% 

Obviously, 

7. A £ °A 

For molecule gaining from corresponding atom recombination: 
I'M = F [ m + (1 - + 7 A*X 


(A = O, N) 


( 1 ) 


(M = 0 2 , N 2 ) 


( 2 ) 


For all other atoms and molecules (surface is assumed to be noncatalytic with respect to them): 


= + (!-*<)*? +*,'*? (3) 

where denotes a convective property such as mass, momentum, or energy. 

Summing over all the species gives, from equations (1), (2), and (3) , the following expression for 
the net flux of momentum or energy: 


NS NS NS NS 

E^ = E^+£( 1 -^+E ( vr < 4 > 


That is, the net flux at the outer edge of the Knudsen layer equals the incident flux, plus the specularly 
reflected flux (incident minus the fraction that sticks) at the wall, plus the diffusely reflected flux 
(those that accommodate to the wall) from the wall. 

Each species is treated separately in the mass balance equations. Therefore, equation (1), (2), or 
(3) is employed, depending on the species being considered, In equation (1), the diffusely reflected 
flux consists of those atoms that are accommodated to the wall minus those that recombine. For the 
molecules in equation (2), the diffusely reflected term is present along with the source term resulting 
from the appropriate atoms recombining on the surface. 

Various Fluxes and Distribution Functions 

The interaction model of equations (1) through (4) is employed to obtain the slip boundary 
equations at the gas-solid interface. Through equation (1) , the net fluxes of species, momentum, and 
energy at the outer edge of the Knudsen layer are equated to the difference between the incident 
and reflected fluxes at the wall. These fluxes are assumed to be constant across the Knudsen layer 
and are obtained from moments of the distribution function. For a convected property <^(V) such 
as mass, momentum, or energy for the ith species, the net flux of that property normal to the wall 
at the outer edge of the Knudsen layer, for example, is 

/ oo roo roo . 

/ / v'y *'(V) /J(V) <i 3 V (5) 

-oo J — OO J— oo 


where Vy is the normal component of the molecular (or total) velocity and /|(V) is the velocity 
distribution function at the edge of the Knudsen layer. 

Similar integrals are obtained for the incident and reflected fluxes by integrations over appropriate 
half-spaces in molecular velocity as shown through the following equations: 

Incident flux: 

F l = f°° f 0 f°° j ^-(V) /l(V) d 3 V (6) 

J — oo J — oo J — oo 
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Specularly reflected flux: 


t roo roo roo ... 

Fj = / vl #(V) f s (V x , -Vy, V z ) d 3 V (7) 

oo JO J— oo 

Diffusely reflected flux; 

/ oo /voo roo . . 

/ / <4 *'CV) A(V) (8) 

-oo JO J— oo 

where /^, is the Maxwellian velocity distribution function evaluated at the wall condition and V x , 
Vy, and V z are the components of the thermal velocity vector V. It may be noted that the normal 

component of the macroscopic (or mass-averaged) velocity VQy at the surface is zero, so that Vy — Vy 
in equations (5) through (8). 

The velocity distribution functions used in integrals contained in relations (5) through (8) are 
those for a nonuniform multicomponent mixture perturbed out of equilibrium: 

f(V) = /( °)*’(V)(1 + **'(VJ] (9) 

where /(0)*(V) is the Maxwellian distribution function for the tth species given as 


and 


/(°)*( V ) 


n»m- 
1 % 


3/2 


(27 rfcT) 3 / 2 



( 10 ) 


®*(V) = 


( 11 ) 


Here k and £ are the dummy indices for three coordinate directions, and the summation with index 
j represents summation over all the species. The summation convention for repeated indices is used. 

The coefficients A B^g, and are functions of the dimensionless velocity: 


W i = M V i (12) 

and are defined as 

4 =[»*> + “a (5 -**?)]»! (13) 

Bit = ‘to (w' k W} - i W?* M ) (14) 

c k i] = 4o )w 'k (13) 

where o^q, fyg, and c) q are constants determined from the variational problem in the first 
approximation for a mixture as given in reference 7, and S^g is the Kronecker delta such that 
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ht = i (* = *) 



6 ki = 0 (fc/£) 

These constants are functions of the collision integrals and are related to the transport properties. 
The simplified form of these constants is provided in appendix A of this paper. More details can 

be found in appendix A of reference 2 or in reference 7. The diffusion vector dj. in equation (11) is 
related to the diffusion velocity of the jth species and is defined after neglecting the external forces 
as (ref. 7): 


d j = JL + ( n l_ W) d]jgp 
k dxfc \ n ) \n p J dxjf. 


(16) 


where «y and my are, respectively, the number density and mass of the jth species; n is the total 
number density, fi = J2j n j\ P is the total mass density, p = J2j n j m j ! an d p is the total pressure, 

p = ^2j pj. A simplified form for d^ is provided in appendix A. 

The total or mass-averaged velocity vg^ (i.e., the kth component) appearing in equation (11) is 
defined as 


NS 


V »k=- 0 'E 


n i m n 


3 = 1 


where vj. is the total velocity (= Vg^ + V^) of the jth species averaged over the distribution 
function, and Vl is the thermal (or peculiar) velocity, also introduced in equation (5). The thermal 
velocity v£ of the jth species averaged over the distribution function is known as the diffusion velocity 



i r 00 roo roo ■ . . 

v k = ^r- / / v l t’M i Y < 17 > 

n j J—OO J—OC3 J—oo 

A simplified expression for the diffusion velocity in terms of transport properties is provided in 
appendix A. 


Balance and Slip Equations 

Based on equations (1) through (4), the balance equations for the *th species for fluxes normal 
to the surface of: (i) the species mass Mj y , ( ii ) the normal component of momentum P^y, (in) the 
tangential component of momentum , and (iv) the energy E^ y are obtained as follows: 

(i) Species mass flux: 

For a recombining atom, 

Ma, = «aM X a + Va - 1a)M I (A = O, N) (18) 

For a molecule gaining from the corresponding atom recombination, 

M My = 0 M M l M + + IaMJ (M = 0 2 , N 2 ) (19) 

For all other atoms and molecules, 


( 20 ) 


M iy = 9 i M} + e i MP 
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(ii) Normal momentum flux: 

NS NS NS . 

E p iv = E < 2 " Wi + E e i p i" = p y < 21 ) 

i—1 i=l i= 1 

(in) Tangential momentmn flux: 

NS NS 

E F iii = E«i p iii < 22 > 

i= 1 i= 1 

(iv) Energy flux: 3 

NS NS 

^ = 53 E iy — 53 E w + 53 e iy M iy 

i—1 i= 1 i 

Diatomic 

molecules 

NS NS 

= E 0 i E P + E + E 

i= 1 i 

Diatomic 

molecules 

+ 53 d t e T M T ( 23 ) 

i 

Diatomic 

molecules 

where is the internal energy of ith species that readily equilibrates with the translational energy 
E 'T under the assumption of local thermodynamic equilibrium. For example, 

rotation __ 

* 

Further, in writing the energy flux balance of equation (23) it has been assumed that there is no 
change in internal energy during specular reflection. 

In obtaining equations (18) through (23) we have used the following relations: 

M i = - M l 4 = 

E 1" = -E^ P% - 0 p( = P? 

% i *11 w ty 

Because it is assumed that the atoms are consumed at the wall by catalytic recombination in 
equation (18), the net mass flux is not equal to zero, M^ y ^ 0. Similarly, the net mass flux is 
not equal to zero, M^ y / 0, in equation (19). However, M^ y = 0, is equal to zero, in equation (20) 
for the atoms and molecules for whom the surface is assumed noncatalytic. 

Substituting equations (5), (6), (7), and (8) with the definition of <^(V) as mass for equations (18) 
through (20), as normal component of momentum for equation (21), as tangential component of 
momentum for equation (22), and as energy for equation (23), respectively, and carrying out the 

3 The energy balance is based on the assumption that the various energies considered readily equilibrate with the 
translational energy. 
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integrations 4 gives equations relating the slip properties to wall properties and gradients at the edge 
of Knudsen layer. All accommodation coefficients 9 j are assumed to be equal to 9. 


Concentration (number density) slip (obtained from mass flux balance): 

/ (Pi,. 


n i _ofi , 1 M iy fZxmj 
nf { ^ 2 nfmiS kT w 



— + 1 

p? ^ 

V l l J 


(24) 


where M^ y and P{y/pf are given by equations (Bl) and (B17), respectively. Alternatively, a more 
simplified expression for Mj y may be obtained from equations (34) and (35). 


Pressure slip (obtained from the flux balance of normal component of momentum) : 

NS 


Ps = <- 


Bh 

'•)/ 


'dv 0x + dv 0z 
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dx 
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L *— 4 J s 
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2 + 


E »* E 4? 4 


3 = 1 


(25) 


with viscosity 


NS 


1 4-^ 

P = zZ n i b iO 
2—1 

Velocity slip (obtained from the flux balance of tangential component of momentum): 

I... X , NS 


4r = 


v 0z = 
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c i 0 
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- S E 
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3= 1 


- 

sj 


i = 1 


dlnT 
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/ NS 
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dVQy 

dz 


\ l JdlnT ( \ 

), + 2 E Pi pr l“ i0 - 2 a>1 ) 


2kT s 


(26) 


(27) 


4 The mass, momentum, and energy fluxes in terms of evaluated integrals over the distribution function are given in 
appendix B. 



Temperature slip; 


(2 kT s f! 2 



NS NS 

Y E iy + (2kT w f/ 2 Y 


i=l 





x 





( 28 ) 


Equations (24) through (28) differ from the corresponding expressions provided in reference 2 be- 
cause of the differences in the interaction model employed at the gas-solid interface and typographical 
errors. 

The constants a^Q, a^, and c^' (also known as the Sonine expansion coefficients) appearing 
in equations (24) through (28) may now be expressed in terms of the transport properties as given in 
appendix A. By using these relations along with the various flux expressions of appendix B and also 

expressing dj. in terms of the gradient of mass fractions dCj/dy as given in appendix A (by neglecting 
diffusion due to pressure gradients and external forces), the following equations are obtained after 
some algebraic simplifications: 5 

Concentration (number density) slip: 

Please note that the following equation is the same as equation (24) and is unchanged: 


Pressure slip: 
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dv 0x dv 0z o dv 0y 
M1 dx + dz 


J S 


2 (2 - 6) / 1 dT ™ 
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2 kT 
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(29) 
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' s 



Velocity Slip: 
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-f 


1/1 dT 


NS 


-i s 


5 \ nkT dx 
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"• E E D a 
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' NS 
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( 30 ) 


5 The approximations made in the expressions for o, i and b t Q are given in appendix A. 
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ft ( dvpz t dv 0y 


VkT V dy dz 


0 


1 / 1 ST 2 ^ 

+ 5laf alE 

}- 1 


- n *E ^ E d vj 

i= 1 i=l 


dz 


C 'S 5 S).)/ 1 '^ 
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Temperature slip: 


ffiu 
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molecules 
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2 kT s m s 
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(32) 


where the mass fraction Q and the mass of a mixture molecule m are defined, respectively, as 




njfnj 

P 


and 



Equations (24), (29), (30), (31), and (32) differ from equations (23) and (25) through (28) of 
reference 3. For dissociated air, the gas-surface interaction model employed in reference 3 appears 


to be inconsistent. Further, the diffusion vector djj. appears to have been incorrectly evaluated in 
reference 3. If one carries out the simplifications in equations (23), (24), and (28) of reference 3 
(which contain the binary assumption) through the evaluation of ^ in terms of mass fraction 
as given by equations (35), (36), and (39) of the same reference, erroneous results are obviously 


NS 

obtained because X) 
3 = 1 



(which is zero by definition) is contained as a factor in several of 


these equations. 



For the first-order recombination at the surface, the following relation between the atom mass 
flux Mj± y and the wall number density may be employed: 


M Ay = - k wA n A m A 


(33) 


where the minus sign indicates that the flux is in the direction opposite to the outward normal, and 
the expression for the rate constant k w A with diffusion and slip is (ref. 2) 6 as follows: 


k wA ~ %4 



(34) 


Here 7^ is the recombination coefficient. 

For a fully catalytic wall (74 = 1), the maximum value of the rate constant k w A is limited by 
the surface temperature. The reaction rate constant for a fully catalytic wall with the gas phase in 
chemical equilibrium is often assumed to be infinity for the sake of simplicity. 

If equation (33) is employed, the net mass flux M^ y appearing in equations (24) and (32) may be 
obtained from 


M Ay ~ ~ k wA n A m A 
M My = k wA n A m A 

M iy = 0 


(A = O, N) 
/ A = O for M = O2 \ 
\A = N for M = N 2 ) 
(For all other atoms and molecules) 


(35a) 

(35b) 

(35c) 


Equation (24), with M iy defined by equations (35), gives the number density ratio nf/nf. However, 
to obtain n™ from this ratio, an expression for n| is required. 

The net mass flux of O and N atoms to the surface, M^ y , is also equal to the rate of consumption 
of these atoms at the wall from surface recombination: 


MAy^'IA^ (A = O, N) 

and the corresponding net mass flux of O2 and N2 molecules will be 


M My = ~lA M Ay 


/ A = O for M = 0 2 \ 
\A = N for M = N 2 / 


(36a) 


(36b) 


For all other species, the net mass flux to the surface may be assumed to be zero. Thus, 

M iy = 0 (36c) 


6 As shown in reference 2, by neglecting slip but keeping diffusion, a slightly different form of equation (34) is obtained: 


kwA — 


2 ^ a I kTw 
2 - nr A V 27rm A 
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Substitution of the net and incident fluxes M^ y and m} from equations (Bl) and (B2), respectively, 


into equation (36) gives the following expression (after thermal diffusion is neglected): 



NS 


m 


dC, 


L ~c i y . 

J Z-^ m q dy 

q = 1 


1 - 




m; 


which may be used to obtain nfj Here ipf is the source term defined as 

« = ia 

A (2 - 74) y/n 2 n s y 

For O 2 and N 2 molecules 

( A = O 

= U-M 


For all other species 


# = 0 


(37) 


H'-H 

(38 a) 

m = o 2 \ 

M = N 2 ) 

(38b) 


(38c) 


Simplifications for Multicomponent Mixture 

Equations (24) and (29) through (32) for multicomponent gas flows can be simplified if one makes 
the following assumptions: 

1. All diffusion coefficients D^j for a multicomponent gas have the same value so that D t j = 

Dam i s the same as the binary diffusion coefficient Dam- 


7 If no assumption is made about Dij, 
species: r 


equation (37) would give the following expression for (dCi/dy) s for all the 

/ \ 1 



m x 

m s 


NS 



\j= l 
L 3& 


ac NS 

q=l 


m dC q 
m q dy 

J s 


mi 


NS 


'Z'&ijCj). 

3 = 1 
3^i 


The source term ipi in this expression, however, may be simplified to yield the following expression for C\ (or n s A ) for 
the recombining atoms only: 


C% = 


Ps 


2VF— lA 7 ^, 

1 A m s 



D S Aj 



NS 


£ 

n— 1 




8 This is a somewhat stronger assumption because is concentration dependent, whereas Ai.? is virtually independent 
of composition. The multicomponent diffusion coefficient D tJ is related to the binary diffusion coefficient through 
the following relations (ref. 7): 

Dij = Kij (-A^j / Mj)Ka 


where quantities Kij are coefficients in a matrix which is the inverse of the matrix with the following coefficients: 

NS _ 


K, 




MqD{q 

q^i 


(* + J) 


— o (* = i) 

One can see, therefore, that by employing Dam (which is the same as Dam) for all the species in a multicomponent 
gas mixture, considerable saving is obtained in computational effort and time without losing the general flavor of 
multicomponent diffusion. This is particularly true if the dissociated air consists predominantly of nitrogen molecules 
and oxygen atoms. The Shuttle entry conditions fall into this category. Simple and multicomponent diffusion gave the 
same results in reference 8. 



2. The ratio of normal momentum flux to the pressure P{ y /Pi is the same for all species and equal 
to that of the mixture; this also implies that the normal shear stress r yy ^ for species i is the same 
as that for the mixture r yy . 

3. The rotational and vibrational states are fully excited so that the internal energy e z - for the 
diatomic molecules in air may be taken as equal to 2 kT/m^ at the edge of Knudsen layer. Implicit 
here is the assumption that excitation of the electronic states of the atoms and molecules does 
not contribute to the internal energy of the particle. This will be true for a relatively cold surface. 

These assumptions retain the major effects of multicomponent fluxes on various slip quantities and 
provide considerable saving in the computational effort required for the analysis of a flow problem 
(refs. 2 and 8). With these simplifications, equations (24) and (29) through (32) yield the following 
results: 

Concentration slip: 

n i =9 (< , 1 M iy l 2 « m i 

nf \ ± 2nfm i ykT w 

where M^ y may be obtained from equations (35) and Py/p s is given by equation (B17) with the 
simplifications mentioned earlier in this section. 

Pressure slip: 
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Velocity slip: 
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Temperature slip: 
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( 43 a) 


or, if the internal energy is frozen during reflection from the surface (see eq. (28) of ref. 5, for 
example), 


Ta _ _ V'' ^iy . 1 ( ?v. a. A V"' / 2fcT s m s _ 3 

r» Z-f m; 2^ + ILV mi mj 

L 1 v y i—l J 


I — 2 — 0 / 1 /C dT 5 M{y \ 

0 \ 2 p dy 4 2^ nmi J 
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Equation (37) may also be simplified to yield an explicit expression for nf: 


where 


n s _ £l < 

M 

\( dCi ) H 

h 1 ( A ^ 1 

/v (- 

acA 1 

Tl — \ 

1 m i 

A dy ) s 

m i \PamJ s. 
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dy)s\ 
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2 - 7 A 2 n s 

^ = -*i 


V 


4>i =o 


(A = O, JV) 

Al = O for M = 0 2 \ 
A = iV for M = N 2 ) 
(For all other species) 


( 43 b) 


(44) 


(45a) 

(45b) 

(45c) 


It is suggested here that the concentration for the major species (for example, nitrogen) be 
obtained by requiring the sum of concentrations of all the species to equal unity. It should be 
mentioned that the mass of the ith species m z - is related to the molecular (or atomic) weight 
through the relation 


m,- k 
Wi~R 


(46) 


where k is the Boltzmann constant and R is the universal gas constant. 
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Concentration Slip Boundary Condition for Fully Catalytic or Noncatalytic Surface 

Equation (44) gives slip values of the concentration ra| for a finite catalytic surface. For a fully 
catalytic (7^4 = 1) surface, one generally assumes complete recombination of atoms at the surface. 
There is a discrepancy in this assumption because the maximum recombination rate is limited by 
the surface temperature as discussed earlier. (See the discussion following equation (34).) Thus, for 
a fully catalytic surface, equation (44) should be employed with 74 = 1. 

For a noncatalytic surface (7^4 = 0), equations (45) give = 0 for all the species. Then 
equation (44) becomes 



which may also be written as 



(47a) 


(47b) 


Summing the above equation over all the species gives 



m dCq\ 
m q dy ) s 


= 0 


(48) 


for a noncatalytic surface. Therefore, equation (47a) is not an appropriate boundary condition. 
However, employing equation (48) with equation (47b) yields 



(49) 


which may be used as the boundary condition for a noncatalytic surface with a multicomponent gas 
mixture. 


Simplifications for Binary Mixture 

At lower altitudes (i.e., when the Reynolds number parameter e 2 of reference 9 is much less than 
1) with the flow in slight nonequilibrium, the derivatives of various quantities with respect to x and 
z and some other higher order terms like dvQ y / dy may be dropped through an order of magnitude 
analysis. (See ref. 10, for example.) Further, the dissociated air may be considered as a binary 
mixture (i.e., consisting of atoms and molecules only, see ref. 9) at these altitudes. If an assumption 
is also made that the internal energy of the molecules remains frozen during reflection from the 
surface, equations (39) through (44) can be simplified further to the following forms. 


Concentration slip: 

n A fT™. 

n\ 2 V T s 

Pressure slip: 
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(51) 
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Velocity slip: 


v 0x = 


^ItHT / [(v^I " v^m) ”A + 


(52) 


u 0* = ^1%-^ ( ^t ) 8 / (V^4 - \/™m) n A + n * 


Temperature slip: 
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In equations (54), is obtained from by using equation (50). 


(54t 


Expression for n^: 
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J_\ SCa 
m M ) dy 



2 7A / W A 

2 - 7A V 27rm A P^AM 


which may also be written as 


n A m A_ r „ _2 - 7A /ssy ns 



s 


(55a) 


(55b) 



For a fully catalytic surface (74 — 1), equation (55b) gives appropriate value for the concentration 
slip. 9 The relevant boundary condition for a noncatalytic surface (7^4 = 0) is 



= 0 


(56) 


as can readily be seen from equation (55b) . 

Comparison With Existing Expressions 

For a noncatalytic surface (74 = 0), pressure and temperature slip equations (51) and (54b) are 
further simplified, with the help of equation (56), to the following equations: 


Pressure slip: 


+ (^Ma + 1~£1 Ku ) 

5x/i 9 \s/Wr dy) s \ ^/W£ M J 


Temperature slip (with frozen internal energy during reflection from the surface): 


(57) 


Tg — 1w 4- y/ft 


2 -9 (K 


29 \k y/2kT 


1 dT\ 
IkT dy). 


m s n s 


1 1_ + JL 

3/2 3/2 1 ^ ^ 3/2 

m M 7 m M J 


(58) 


Equations (52), (53), (57), and (58) ( with the noncatalytic surface assumption in equations (57) and 
(58)) can be rewritten as 


Velocity slip: 


* " Aiai n s Vm s (vs 5 «r) 8 


(59) 


Pressure slip: 


Temperature slip: 


v s = ( JL 

1 I sV / Ps WfcT dy J t 


p s = p w + Bib 1 


(V) 



— 

s x • kT dy 


Ts -T W + Cici( 1 ^ ) g (fc VfcT ty/vm, 


(60) 


(61) 


(62) 


9 Sometimes in the literature Ca is prescribed as zero for a fully catalytic surface. Strictly speaking, this is true only 
when the Reynolds number parameter e 2 (ref. 9) is approximately zero (close to the chemical equilibrium condition at 
low altitudes) and the surface is highly cooled. This can easily be seen by nondimensionalizing equation (55b) in a way 
similar to equation (2.8d) of reference 9. Thus, the recombination rate coefficient 7 a and density (as measured through 
£ 2 ) both control the recombination rate and not 7 a alone. 
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where 
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nsVrns 
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0 
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0 

2-0 
0 
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0 0 

The concentration slip condition consistent with equations (59), (60), (61), and (62) is 



(63) 


The range of values of /Ip Up and C\ is as follows for a mixture of oxygen atoms and molecules and 
7s = 1.4 (ref. 9): 

1.0039 < Ai < 1.0186 

0.9507 < Bi < 0.9507 ( J + ./? Km \ j ( K A K M 

~ y V 4 m A + V 2 m M J / + ^ 

0.9056 < Ci < 0.9507 

The minimum values for Ai and C\ occur at = O.bp/m 4, whereas the maximum value for -Bj 

occurs at this value of n^» 

Expressions (59), (61), and (62) reduce to those obtained in reference 9 10 if one assumes the 
values for Tj, -Bj, and C\ to be unity and a noncatalytic surface boimdary condition as given by 
equation (63). The concentration slip boundary condition provided in reference 9 is 




VL+ 


2 - 1A 

2 1A 




kT s 


-P S AM 


(9Ca\ 

V dn* J , 


(64) 


which is the same as equation (55b) obtained herein except for the additive term C^. 11 It is 
obvious that the temperature slip boundary condition of equation (2.8c) obtained in reference 9 is 


10 See appendix C for the dimensional form of the slip boundary conditions given in reference 9. 

11 It appears that equation (2.8d) of reference 9 for the concentration slip has been formulated for the chemical 

equilibrium wall condition. Thus, when surface slip goes to zero, one obtains from equation (2.8d) the equilibrium value 
for the concentration, i.e., which would be zero for the oxygen atoms for surface temperatures of 2000 K 

or less. For a highly cooled wall, the equilibrium catalytic wall condition can be approximated by a fully catalytic wall 
condition where the species concentrations are set equal to the cold gas equilibrium value. 
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valid, strictly speaking, for a noncatalytic surface only and is not consistent with the concentration 
slip boundary condition of equation (2.8d) (reproduced here as equation (64)) obtained for a finite 
catalytic wall. As a matter of fact, the velocity and temperature slip boundary conditions of 
reference 9 are similar to those of reference 10, 12 where these slip conditions are provided for a 
perfect gas (or single-species gas). 

The inconsistencies in the boundary conditions used in references 11 and 12 are similar to those 
of reference 10, namely, the pressure and temperature slip values are in error by the ratio of 
dimensionless thermal conductivity to the viscosity K/p. References 11 and 12 also employ the 
single-species formulation for obtaining the velocity, temperature, and pressure slip values, and the 
concentration slip expression in these references is the one obtained for a binary mixture (ref. 9). 
The coefficients appearing in the slip equations of reference 12 may be obtained by multiplying op, 
&1, and c\ (coefficients of ref. 9) given here by 16/57T. The factor 16/5 is sometimes replaced by 7r. 

Simplifications for Single-Species Gas 

Equations (39) through (42) and equation (43b) may be simplified for a single-species gas to the 
expressions obtained in reference 6. For single species, these simplified equations may be written as 
(with 7^4 = 0 and M^ y = 0 for a single-species gas). 


Density slip: 
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mn s 

jTw 

[i + 5 4 1* A « | 

( dv 0x 

, dv 0z 

J v 0 y\ 1 
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mfi w ^ 

/ T s 

L ^ 24 V 2 ^RT S ' 

\ dx 

H dz 

J dy ) s. 


Pressure slip: 



Temperature slip: 13 


Ts_ 

Tw 



75-7T 2 — 9 
128 0 



(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


12 There appears to be some error with the form of equations (2.7c) and (2.7d) given in reference 10 if one employs 
the definition from equation (2.4b) of the same reference of dimensionless heat- transfer rate q in these equations. (The 
Reynolds number parameter e 1 is missing in eqs. (2.2b) and (2.4b) of the same reference for dimensionless shear stress 
and heat transfer, respectively. These have been corrected in ref. 9.) The resulting pressure and temperature slip 
expressions contain the dimensionless viscosity coefficient fi rather than the dimensionless thermal conductivity K. 
References 11 and 12 contain same error. See equation (2.8c) of reference 9 for the “correct” form. 

13 Since the temperature slip expression contained in reference 6 is for a gas consisting of perfectly elastic spherical 
molecules, the gas does not possess any internal degree of freedom. Therefore, we have obtained the temperature slip 
equation from equation (43b). 
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where we have used the following relations between the coefficients of viscosity and thermal 
conductivity and the mean free path for perfectly elastic spherical molecules (ref. 13): 


At = jg p\s/2ttRT 

v 15 k 
K = — =p 
4 m 


(70a) 


(70b) 


The expression for /t given earlier, following equation (25), for a single-species mixture is modified to 

(70c) 


l =fm 
At = -nkTb 0 


Now, if an assumption is made that the slip values are small , equations for the density, pressure, and 
temperature slip are simplified to 


Density slip: 


Ps_ 

Pw 



1 _ A f* hi ( 

24 V 2 ^RT^\ 


dvpx + dvQz _ 2 dv 0y 


dx 


dz 


dy 


)J 


(71) 


Pressure slip: 
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(73a) 


(73b) 


In obtaining equation (73b) from equation (73a), we have used Py/p s ss 1 in the second-order terms 
(i.e., 2d and 3d terms on right side of eq. (73a)). Equations (67), (68), (72), and (73b) are the slip 
equations given in references 6 and 13. 


No-Slip Boundary Condition for Species Concentration 
Multicomponent Mixture 

The no-slip boundary condition for the species concentration without any assumptions may be 
obtained from equations (44) and (45). In the absence of slip, the thickness of the Knudsen layer 
shrinks to zero and the values at the top of the Knudsen layer become the wall values (see fig. 1): 
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with 
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(For all other species) 
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where we have neglected the higher order shear (i.e., Py/p w ~ 1)- 

It is suggested here, again, that the concentration for the major species (for example, nitrogen) 
be obtained by requiring the sum of concentrations of all the species to equal unity. 

For the recombining O and N atoms, equation (74) may further be simplified to 


n A = 17— I nV AM 
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NS = 

m dUq 
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(76) 


The recombination rate constant k wA in equation (76) has been defined as (ref. 2) 

27 A I 


^wA 


w 


2 - 1A V 2nm A 


without slip and higher order shear. 

For a noncatalytic wall (with 74 = 0), equation (49) gives 

©)„=• 

for all the species of a multicomponent mixture. 

Binary Mixture 

For a two-species mixture of atoms and molecules, equation (55b) gives 


r> w _ _ 

A Ka 
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D kU 


9Ca \ 

dy A 


for a surface with finite catalyticity and 


( 7 ).=" 


(77) 


(78) 


(79) 


(80) 


for a noncatalytic surface. Equations (79) and (80) are used in reference 4. 

Appendixes D and E give the slip and no-slip boundary conditions (presented in the text earlier), in 
the dimensionless form for the body-fitted and spherical polar coordinates, respectively, for a planar 
flow. The various integrals employed for evaluation of the net, incident, and specularly reflected 
fluxes defined by equations (5), (6), and (7), respectively, are given in appendix F. 
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Discussion 

The present analysis provides a consistent formulation for the slip equations for multicomponent 
and binary mixtures and single-species gas reacting catalytically on the surface. The diff erences 
between the slip expressions obtained by various researchers have been reconciled, and the implication 
of various assumptions (some of them inconsistent) contained in those expressions is discussed. The 
slip equations have also been obtained in body-fitted and spherical polar coordinates in a form which 
can readily be employed in the flow-field calculations. 

Usually, the equation for pressure slip is not required as a boundary condition but is needed to 
obtain the surface pressure. The temperature slip equation given herein is for a constant surface 
temperature, which is provided as a boundary condition. For an adiabatic surface, however, the slip 
temperature T s may be obtained by equating the wall heat-transfer rate to zero; that is 

^=( K f-E-'A+«^) s = o 

where the expression for mass flux J ? - is provided in appendix B and the higher order terms have 
been dropped. The temperature slip equation is now required to obtain the wall temperature T w . 

An expression has also been obtained for the finite-rate species-concentration boundary condition 
for a multicomponent gas mixture without surface slip. This boundary condition in the literature 
(ref. 4) has generally been specified by assuming the dissociated air to be a binary mixture. However, 
the binary mixture formulation does not provide boundary conditions for the recombined molecules 
(on the surface) such as O 2 and N 2 in a multicomponent mixture; it gives boundary conditions only 
for the recombining atoms. 

To assess the importance of various terms in the general boundary conditions suggested herein, 
a study is being undertaken in which the flow properties and boundary conditions would be varied 
systematically. The boundary equations form a simultaneous set, which is being coupled with a 
flow-field calculation procedure in the stagnation region. These equations would finally be coupled 
with the viscous shock-layer code developed by Moss (ref. 14) for the detailed flow-field calculations. 
The boundary equations obtained in the present work should provide a more realistic set of boundary 
conditions for a multicomponent mixture for low Reynolds number slip flows as well as no-slip flows. 

Concluding Remarks 

The boundary slip expressions obtained in this paper are closed-form solutions of the mass, 
momentum, and energy flux equations by using the Chapman-Enskog velocity distribution function. 
This function represents a solution of the Boltzmann equation in the Navier-Stokes approximation. 
The analysis, obtained for nonequilibrium multicomponent airflow, includes the finite-rate surface 
catalytic recombination and changes in the internal energy during reflection from the surface. 
Expressions for the various slip quantities have been obtained in a form which can readily be 
employed in flow-field computations. A consistent set of equations is provided for multicomponent 
and binary mixtures and single-species gas. An expression is also provided for the finite-rate species- 
concentration boundary condition for a multicomponent mixture in the absence of slip. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
May 17, 1985 
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Appendix A 


Constants Appearing in Distribution Function and Diffusion Velocity in Terms of 
Transport Properties 

General Expressions for Sonine Expansion Coefficients 

(j\ 

The Sonine expansion coefficients a z -g, a ?1 , b^Q, and c^' appearing in the general velocity- 
distribution function are found by a variational technique in which they are solutions to sets of 
simultaneous equations. References 2 and 7 provide the coefficients in terms of solutions to these 

(s l) 

sets of equations. These solutions are expressed in terms of the collision integrals 0 Y ' . Reference 7 
also provides the transport properties in terms of the Sonine expansion coefficients. Thus, in place 
of evaluating these coefficients in terms of the collision integrals, they may be expressed in terms of 
the transport properties. The various relations are 
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polynomial coefficients . 


Now cj^ = 0 (see ref. 7), then 
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(A4) 


where Dj(£), K(£), /^(£), and jD^-(£) are the multicomponent thermal diffusion coefficient, thermal 
conductivity, coefficient of viscosity, and the multicomponent diffusion coefficient, respectively. The 
argument £ is the number of terms used in the Sonine expansion. Except for Dj, letting £ = 1 

gives quite good results for K, p, and D^. When £ = 1, however, coefficients Dj vanish. Hence, 
in order to get the coefficient of thermal diffusion, it is necessary to take at least two terms in the 
Sonine expansion (i.e., £ = 2). If the argument does not appear with a coefficient (except for g), it 
is considered to have one term in the Sonine expansion. 

The diffusion velocity for ith species in terms of the transport coefficients is obtained as (see 
ref. 7) 



where the diffusion vector d J k for the jth species is defined after neglecting the external forces as (see 
eq. (16)): 
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Equation (A6) may be simplified to obtain 
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Simplified Expressions for Coefficients a tl and 

Since an and 6 z0 cannot be obtained directly from equations (A2) and (A3), we introduce the 
following assumptions for the mixture thermal conductivity and viscosity, respectively, 


NS 

K *T, % K i 

i= 1 

(A8) 

NS 

** « E w Ui 

i= 1 

(A9) 

Equations (A8) and (A9) are approximate forms for the more exact formula of Wilke. (See ref. 15.) 
These equations imply that 

n j 

E * 1 

3 = 1 

(A10) 

in Wilke’s relation. This relation is approximately true for air. 

Thus, with the help of equations (A2), (A3), (A8), and (A9), 
approximate expressions for an and respectively, 

we may obtain the following 

4 rij Ka I rriA 
k^kT 

(All) 

i 2 rii 

"A»“i irs^ 

(A12) 


For accurate evaluations of an and 6 z q, cumbersome expressions of the type of equation (7.4-49) of 
reference 7 are required to be solved. 
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Appendix B 


Evaluated Expressions for Fluxes of Mass, Momentum, and Energy From Integrals 
Over Velocity Distribution Function 

The fluxes of mass, momentum, and energy (i.e., of a property <£ Z (V)) are given in terms of 
integrals over the velocity distribution function /*(V) in the section “Analysis” of the main text. 
Here we give evaluated forms for these fluxes obtained from integration over the distribution functions. 
The various integrals needed in these evaluations are provided in appendix F. 

(i) Expressions for mass flux of species: 


Net 


Incident 


^ik — 
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m,n| ^y/2 ^ dlnr _ 


"2 V 
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a i0 + n E c iO d k 
3- 1 


Specularly reflected 
Diffusely reflected 


/2fcr s y/ 2 
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'9v 0x 
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V m i > 
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(k = x, y, z) 


-M:. 


(Bl) 


(B2) 

(B3) 

(B4) 


It may be mentioned that M^ y as obtained here is similar to the mass flux Jj expression of reference 16. 
(ii) Correlation with Blottner’s expression (ref. 16) for mass flux: 

If the thermal diffusion term is neglected in equation (Bl) and substitutions are made for 
and di from equations (A4) and (A7), equation (Bl) may be simplified to obtain 14 


— Afjy — 


pWi 
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NS 
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D, 


dC; = 


13 


dy 


CjW 


NS 

E 

Q = 1 


1 dCq 

Wq~dy~ 


(B5) 


where is the molecular (or atomic) weight of the *th species and is related to the mass of the ith 
species m ^ through the relation 


m,- k 

Wi~J 

Here k is the Boltzmann constant and R is the "universal gas constant. Equation (B5) may also be 
written as 



\y. NS 

Ji=pU y 

w h 


D, 


w NS 

ij ~ yf. E D iq C< l 
3 qfr 


dOj 

dy 


14 It may be noted here that Du = 0. See reference 7 for details. 


(B7a) 
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or, with the introduction of the multicomponent Lewis number Np e defined as 


_ P c pDjj pNp r Djj 

iy Le,ij ~ ~ K — - 
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Equation (B7a) may be rewritten as 
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In general, 15 


^ Le,ij ^Le,ij (^t , , Mj , 

If the binary Lewis numbers ^i e ^j are assumed to be the same for all the species, then 

^ Le,ij ~ N Le,ij = ^ Le,ij(Q’^v M j) 
and equation (B8) for A fyy becomes 

\ NS 


M*ij = N Le,i ~ 


W 


V 3 J q^i 


and the mass flux (eq. (B7e)) now becomes 
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(B7b) 
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(B 8 ) 
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(BIO) 


15 


See equations (7) and (8) of reference 16. 



Further, if the binary Lewis numbers ^Le,ij are constant for all the species, the term 


NS = * dC, 
A b*- — J - 
*3 dy 


E 

3& 


is zero in equation (BIO). Therefore, 
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} (S)" 


(Blla) 


(Bllb) 


Thus, equation (BIO) gives the mass flux due to concentration gradients only with A b*- and Ni e j 
defined by equations (B9) and (Bllb), respectively. 

To utilize equations (B9) and (Bllb), N^ e ^ is still to be evaluated through complex matrix 
inversions (ref. 16). If we now make a little stronger assumption such that iV£ g - = N^am = N^ e 
(same for all the species), we obtain a much simpler expression for A6| -, and equation (Bllb) becomes 


or 
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(B12a) 


(B12b) 


Similarly the expression for A 6* • may be obtained from equation (B9) as 

1 IJ 


A % = N Le4 - N Le 




(B13) 


It may be noted that for — Q, the term Ni e j in equation (BIO) vanishes and equa- 
tion (Bllb) or (B12b) is not required. Further, employing a constant value for the Lewis number 
N Le P^PtBam . j , d oes no t imply that has a constant value. 

It should be pointed out here that for the case when J^Le,ij the same for all the species, we 
have used the condition 


NS 


E Si 4 

.3# 


dCj_ 

dy 


= 0 


(and not just A6*^. = 0 as employed in ref. 16) to obtain equation (Bllb) for iVjr^. The present 
formulation appears to provide a more consistent expression for N^ e ^. 
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Finally, equations (BIO), (B12b), and (B13) for a binary (two-species) mixture give 


which is Fick’s law for mass flux 


j. = -JL~ Nt && 

1 Np r Le dy 


N Le,i = N Le 


= 0 


(in) Expressions for normal momentum flux: 
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(iv) Expressions for tangential momentum flux: 


Net 


P ix ~ 


~Pi 


P iz = ~P S i 


o| 

•<s> 

*ol 

( dv 0x 
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p j- = _p! = _L p? 


ldlnT 
+ 2 dz 
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ix zz u 


(v) Expressions for energy flux: 
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Appendix C 

Slip Boundary Conditions Given by Davis for Binary Mixture 

In this appendix, the dimensional form of the slip boundary conditions obtained by Davis (refs. 9 
and 10) is provided for the purpose of comparison with the slip expressions obtained in the main 
text for a binary mixture. 

The dimensionless form of the surface slip conditions provided in references 9 and 10 are 
Concentration slip (ref. 9): 16 




C!,e + * 2 


2-7^4 {F D ) 8 
2 W %4 Ps 



(Cl) 


Velocity slip (ref. 9): 


Pressure slip (ref. 10) : 17 


Temperature slip (ref. 9): 


u s 


_ 2 Ps Fps_ f du \ 
~ £ ai Ps \ Ps \dn) s 


__ i £.2 l K s Ps(dT\ 

- Pw+£h T s yp s \dn) t 


2 K s lp s (dT\ 
T w + e Cl _^/-( — 


m 


(C3) 


(C4) 


With the help of the nondimensionalizing quantities given in references 9 and 10, the following 
dimensional forms of equations (Cl) through (C4) are obtained: 


Concentration slip: 




(C5) 

Velocity slip: 

Ps ( 1 

Ua ~ ai y/m\nVm 

(C6) 

Pressure slip: 




Temperature slip: 



16 Here the dimensionless quantities are denoted with a bar over the symbol. See appendix D for details of 
nondimensionalization and the coordinate system. 

17 For the reasons explained in the main text, we have employed K in place of ji in this expression. 
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where n* is the coordinate normal to the surface (see fig. Dl), and n and m are the number density 
and mass of a mixture molecule, respectively. The constants aj, b\, and c\ are defined as 


«1 = 


h = Tz 


c i 



2-0 2 
— — = 1.2304- 
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15 Hr 2 — 9 

16 V 2 0 - 1 * 1750 
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" 128 V 2’ 
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= 2.3071 
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2-0 

0 

2-0 


0 
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Appendix D 

Slip Conditions in Body-Oriented Coordinate System 

In this appendix, we obtain the slip conditions for an axisymraetric body in the body-oriented 
coordinate (s*,n*) system. The coordinate configuration for such a system is given in figure Dl. The 
equations of this appendix contain simplifications similar to those given in the section “Simplifications 
for Multicomponent Mixture” in the main text. By employing the metric coefficients (ref. 17) of 
hi = 1 + n*K and «= 1, equations (39) through (45) give the following slip equations in the 
body-oriented coordinate system. 



Figure Dl. Coordinate configuration for body that is symmetric about its axis. 
Concentration slip: 
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where (from eqs. (35)) we define Mj y as 
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Velocity slip: 
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If the internal energy (comprising rotational and vibrational energies in this case) is considered frozen 
during reflection from the surface, the following form of the temperature slip equation is obtained 
(see eq. (28) of ref. 5): 



Equation for n|: 



where 
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In equations (Dl), (D5), (D6), and (D8a), the ratio 
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and 6g is related to viscosity p through the relation 
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(D8b) 

(D8c) 
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which has been obtained from equation (A3) by assuming to be the same (6g) for all the species. 

Equations (Dl) through (D8) may now be nondimensionalized by using the following relations 
(refs. 9 and 14): 
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where rjy is the body nose radius of curvature, T re f = r/Jb/cp )0 o, and >u re f = p{U%o/cp, oo)- If 
the nondimensional quantities as defined here are used, the following equations are obtained from 
equations (Dl) through (D8): 

Concentration slip: 
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where 
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Temperature slip: 
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or, with frozen internal energy during reflection from surface, 
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where 
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For a noncatalytic surface (74 = kyyA ~ 0), equations (D18) give 'ip? = 0 for all species, 
case, equation (D17) becomes 


(St) 

Equation (D20), similar to equation (47b), gives 
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which may be employed as the boundary condition for a noncatalytic surface with a multicomponent 
gas mixture (with the binary assumption for diffusion coefficients). For a fully catalytic surface, 
employing 74 = 1 and using equations (D18) in equation (D17) would yield the appropriate 
concentration slip. 

If the multicomponent diffusion coefficients are retained without the binary assumption, the 
underlined terms in equation (D13) are to be replaced by 

® NLeJi /££>„ ^ Wac, 

2^ Nr e \ dn 3 Z^Wn dn 
and the underlined terms in equation (D14) are to be replaced by 
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where the multicomponent Lewis number is defined as 
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Simplifications for a Binary Mixture 

When all the species in a gas mixture can be considered as atoms and molecules only (ref. 9) , 
equations (Dll) through (D18) may be further simplified. With the assumption of Py/p s 1, 
equations (Dll) through (D18) for a binary mixture yield 


Concentration slip: 
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where we have used the following equation of state 
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and neglected the higher order shear terms as well as diffusion terms. 



Velocity slip: 
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where we have again employed the equation of state (eq. (D24)) and neglected higher order shear, 
conduction, and diffusion terms. In obtaining equation (D25), we have also used the following 
approximation: 
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For the equation to be more exact, one may keep this factor in equation (D25), 
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where 
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and | < Xl < 1 for 0 < C 4 < 1 and diatomic molecules (W^ = 2 IV 4 ). In obtaining equation (D26), 
we have employed the equation of state (eq. (D24)) and = 2W^. 

With frozen internal energy during reflection from the surface, the temperature slip equation 
becomes 
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For diatomic molecules {Wj^f = 2W^), X2 ma Y t* e taken as unity for 0 < C 4 < 1 . Equation (D27) 
employs Wj^ = 2W^. 

Equation for C%: 


Ca = 


2 - JA _ £ 
2 


(NlA (p\ [j[w (dCA) 
\ Npr ) s \ pJs \ T s T Te{ V dn ) t 


(D28) 


38 



Equation (D28), for a noncatalytic surface with 7^4 = k w A = 0, gives 



whereas using 7,4 = 1 in equation (D28) and in the expression for k Wj 4, 


(D29) 
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would yield the appropriate concentration slip for a fully catalytic surface. 

With the following (somewhat inconsistent) assumptions, equations (D23), (D25), (D27), and 
(D28) may be simplified to those obtained in reference 9 and employed in reference 12: 

(i) In pressure slip equation (D23), N Pr — 1 is employed along with approximations 
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and 



« Vs 


(it) In velocity slip equation (D25), the following assumption is made for the mixture molecular 
weights 

^00 t 

W s 

(in) In temperature slip equation (D27), an assumption is made that the surface is noncatalytic 
so that k w A — 0. In addition, Prandtl number Np r is taken as unity, and it is assumed that 


2 \7-lJ s 8 

(■ iv ) No assumption is required in the concentration slip equation (D28). But it may not be 
appropriate to employ (ref. 12) this equation for all the species of a multicomponent mixture. This 
equation is applicable to concentration slip of recombining atoms only in a binary (two-species) 
mixture (ref. 9) of atoms and molecules. 

No equation has been obtained in reference 9 or 12 to correspond to equation (D22) to obtain 
wall values of the species concentration (C^) from the values at the edge of the Knudsen layer (C^). 


Slip Expressions For Single-Species Gas 


For a single-species gas (74 = 0), the following slip expressions are obtained: 
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Velocity slip: 
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Temperature slip (for a gas consisting of molecules only (i.e., diatomic perfect gas)): 
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where we have used xi — j and assumed that the diatomic molecules possess full rotational as well 
as vibrational degrees of freedom at the edge of the Knudsen layer, implying is m 9/7. Generally 
this would be true due to the high temperatures at the Knudsen layer edge. However, for perfect air, 
7oo is generally taken as 7/5 with the excitation of only the rotational degree of freedom. For such 
an analysis, therefore, 7§ would also be equal to 7/5. In this case (with only the rotational degree of 
freedom excited), equation (D34a) would be modified to 
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With frozen internal energy during reflection from the surface or for a monoatomic gas, we obtain 
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where 75 = 7/5 for perfect air, and 7$ = 5/3 for a monoatomic gas. Equations (D32), (D33), and 
(D35) are the ones employed in references 9 and 12 with Np r = 1. 


No-Slip Species Concentration Boundary Condition 

Multicomponent mixture. The no-slip boundary condition may be obtained from equations (D17) 
and (D18). In the absence of slip, the thickness of the Knudsen layer shrinks to almost zero, and the 
values at the top of the Knudsen layer become the wall values (see fig. 1): 
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and we have neglected the higher order shear (i.e., Py/p w & 1). 

For the recombining O and N atoms, equation (D17) may also be written as 
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The recombination rate constant k wA in equations (D37a) and (D38) has been defined as (ref. 2) 


% — 1 J W cqT w 

wA 2 - Moo V 27T700 W A Too 


without slip and higher order shear. 

For a noncatalytic wall (with 7^4 = k wA = 0), equation (D21) gives 



(D39) 


(D40) 


for all the species of a multicomponent mixture with the binary assumption about the diffusion 
coefficients. 

Binary mixture. For a two-species mixture of atoms and molecules, equation (D28) gives 


°A 




(D41) 


which may also be obtained from equation (D38) for a surface with finite catalyticity. In equa- 
tion (D41), k wA is again obtained from equation (D39). For a noncatalytic surface with k wA — 0, 
equation (D41) gives 



(D42) 


Equations (D41) and (D42) are obtained in reference 4. However, in reference d, equation (D30) is 
employed in place of equation (D39) for obtaining the reaction rate constant k wA . As obtained in 
reference 2, equation (D30) is more appropriate when slip and higher order shear are included in the 
flow-field analysis (e.g., under rarefied or low density conditions). 

If k wA is substituted from equation (D39) into equation (D41), we may also obtain 
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where we have used the following equation of state: 
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Equation (D43) compares with the corresponding equation of reference 9, if one keeps in mind 
that the diffusion coefficient Fjj (of ref. 9) in the absence of slip is related to the present variables 
by the relation 
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It may be seen from equation (D41) or (D43) that the gradient {dC A / dn) w is governed by the 
ratio k wA /e ^ or 1 A /^ ■ For surfaces which are almost noncatalytic, this ratio would be of the order 
of 1 for large values of the Reynolds number parameter ( 1 /e^). This would imply that (dC A /dn)w 
is not necessarily zero for such a surface if ~ 0. In addition, equations (D41) and (D43) also imply 
that regardless of the surface catalytic activity, the higher density conditions (smaller values of e^) 
would promote larger values of the concentration gradient at the surface. 
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Appendix E 

Slip Conditions in Spherical Coordinate System 

For analyzing the flow in the stagnation region of a body, it is convenient to use the spherical 
polar coordinates (r, (/>) for the two-dimensional flow, (See fig. El.) The following relations exist 
between the spherical polar and the body-oriented (s*,n*) coordinates over the spherical portion of 
a body: 

r sin <f> = R* + n* cos (3 

r = r N + n * — r iv(l + n* k) 

1 

K — — 

r N 

rj^<f> = s* 

/3 = 90 — (j> 

or, in the nondimensional form, 

fsin (j> — R + ncoscf) 

r = ?N + n = l + n 

k = 1 

(f> = s 

df — dn 
d4> = ds 

where we have used nose radius rj y to nondimensionalize all the distances. (See appendix D.) 





Figure El. Spherical coordinate configuration. 
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By using the relations given by equation (E2) in equations (Dll) through (D19), we can obtain 
the slip equations in the spherical polar coordinates for a multicomponent mixture. Once again, 
these equations employ the simplifications given in the section “Simplifications for Multicomponent 
Mixture” in the main test. 
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Pressure slip: 
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Velocity slip: 
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Temperature slip: 
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or, with frozen internal energy during reflection from surface, 
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In obtaining these various equations, the following form of the equation of state has been employed: 
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Similar to Appendix D, equations (E3) through (Ell) may further be simplified for a binary 
mixture and for a single-species gas. 



Appendix F 

Integrals Required for Obtaining Various Fluxes 

In order to obtain the normal luxes of mass, momentum, and energy, it is necessary to evaluate 
various integrals over the velocity space of the distribution function. These integrals involve terms 
that are various velocity moments of the distribution function. The integrals are provided for the 
net, incident, and specularly reflected fluxes. These integrals are basically the same as those provided 
in reference 2 except for some corrections and additions. 

Integrals Over Entire Velocity Space (for Net Quantities) 
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Integrals Over Lower Half of Velocity Space (for Incident Quantities) 
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Integrals Over Upper Half of Velocity Space (for Specularly Reflected Quantities) 
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j°° J™ j°° Wy{-Wy)W A exp(-W 2 ) d z W = -||^r 3/2 

The term [—W y ) appearing in the integrals emphasizes that the sign of the thermal velocity 
component normal to the surface V y (and consequently that of W y ) changes in the distribu- 
tion function /j upon reflection from the surface. This implies that the distribution function 
fhv x ,Vy,V z ) = f l s {Vx,-Vy,V z ). 
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